Journal of Phycology

John Raven pioneered the field of algae ecophysiology, advancing our understanding of cellular resource economics, carbon acquisition, and energy allocation. His work laid the foundation for investigating integrative physiology, linking growth-survival trade-offs across diverse environments. The sea ice habitat provides an excellent framework to continue the research John championed. With steep temperature-salinity gradients, algae survival requires a shift in physiology that we are only beginning to understand. We developed two small scale, reproducible icecosms to investigate physiological changes associated with incorporation into sea ice and survival potential post-melt. Fragilariopsis cylindrus and Nitzschia frigida, known for their association with the ice environment, and Porosira glacialis, known for its association with the ice edge, were used to mechanistically link physical properties with algal physiology and post-melt survival. We observe incorporation into the ice of F. cylindrus and N. frigida alongside vertical photophysiological profiles of F. cylindrus revealing inhospitable conditions in the top compared to the bottom layers of ice. N. frigida and P. glacialis remain viable within the ice and retain the capacity to seed populations following melt. Our results establish icecosms as experimental framework to investigate ecophysiological responses of sea ice algae and provide a foundation toward ecological and evolutionary questions.

Photosymbiosis has evolved multiple times independently in ciliates. However, these associations can be antagonized by shifts in environmental parameters that impose stress on the host, necessitating the evolution of mechanisms to contend with this stress and to control the symbiont population. To investigate whether convergent strategies have evolved among algae-bearing ciliates in the class Oligohymenophorea, we imposed light stress on three host species that represent at least two independent evolutionary origins of photosymbiosis and measured their cellular responses. Under high light, all three species experienced an initial drop in host cell density which recovered to levels commensurate with those under low-light conditions as they decreased their symbiont loads. We then performed a comparative transcriptomic study to investigate whether a core set of genes exists that is involved in this response. Thirty-one gene families possess differentially expressed transcripts across all three species that included the upregulation C1 and S28 class peptidases, genes involved in ROS mitigation, and a gene with potential involvement in mitochondrial remodeling associated with changes in algal symbiont load. We additionally found downregulation in Dicer, which could mitigate the processing of algal transcripts by the hosts RNAi machinery that are freed upon algal digestion, and downregulation of motor proteins that may reflect changes in the hosts swimming behaviors and transport of intracellular vesicles in response to light. The 31 gene families are present and widespread in non-symbiotic oligohymenophoreans, illustrating that a pre-existing genetic toolkit exists in this clade that helps explain how it is predisposed to evolving photosymbioses.

Algae require specific acclimation strategies to cope with spectral variability in shallow marine habitats. We investigated how the siphonous green alga Derbesia alters its photosynthetic and metabolic processes under white (WL), blue (BL), green (GL), red (RL), and far-red light (FL) by conducting photobiological and transcriptomic sampling over a 10-day period. Our results show two contrasting photoacclimation strategies: BL and GL promoted metabolic activity associated with growth, whereas FL and RL induced a low-light-like survival strategy characterized by reduced growth and suppression of the core metabolism. Photosynthetic acclimation across all conditions primarily occurs within the light dependent reactions. BL and GL promoted early acclimation marked by the immediate activation of light-harvesting complexes (LHCs) and a key transcriptional regulator MYB, and showed better acclimation marked by the sustained activation of ATPases, ATP transporters, and hormone-signaling components. BL induced a distinct transcriptional shift during the transition to prolonged exposure, including enhanced cyclic electron transport, and key regulators of protein synthesis, DNA replication, and transcriptional regulation. In contrast, FL, and to a lesser extent RL, triggered responses resembling low light acclimation with constrained growth, characterized by inefficient energy utilization, enlarged antenna systems, chloroplast proliferation with aggregations, and reduced growth rates. This study suggests high accumulation of core photopigments and reduction in chlorophyll a/b is an acclimatory response to FL, and consistently higher activation of core metabolic processes under WL likely indicates the evolutionary adaptation of Derbesia to shallow coastal environments where broad-spectrum light predominates. Additionally, our newly sequenced draft genome of the Derbesia strain for this study could serve as a genomic resource for future molecular photobiology research in Bryopsidales algae.

Labyrinthulomycetes are a class of fungus-like heterotrophic protists from the Stramenopiles lineage, recognized for their ecological role as decomposers and contributors to nutrient cycling. They colonize various substrates, from seaweed to terrestrial environments, utilizing ectoplasmic networks for nutrient absorption. This study characterized a novel Labyrinthula strain associated with the marine diatom Biddulphia. Phylogenetic analysis of the full-length 18S rRNA gene positioned this strain as a new species, Labyrinthula merlionensis sp. nov. Scanning electron and light microscopy observations revealed bi-flagellated zoospores and spindle-shaped vegetative cells with ectoplasmic networks. Time-series observations of the interactions between L. merlionensis and Biddulphia were categorised into different phases: establishment, infection, and aggregation. Scanning electron and confocal microscopy observations during the infection phase established the use of ectoplasmic nets to target the marginal ridge regions between diatoms, and the detection of labyrinthulid cells within diatom frustules. These findings enhance the understanding of the diversity, morphology, and ecological roles of Labyrinthulomycetes, particularly their intra- and extra-cellular interactions with diatom hosts.

Diatom mating activity contributes to their enormous phenotypic and genetic diversity, yet little is known about patterns in diatom reproductive compatibility across genetically diverse strains, nor the effects on offspring phenotypes that may confer adaptive evolution, niche partitioning, or trait improvement. Here a panel of 38 Arctic Cylindrotheca sp. isolates were crossed pairwise to detect mating compatibility. Positive mating patterns were identified in multiple clades, including amongst crosses of different parental rbcL genotypes. F1 isolated from three different crosses presented phenotypic variation in growth rate, plastid traits, and associated photo-physiological responses to blue and green actinic light. Offspring gliding speed and behaviour also varied, providing insights into complex motility traits that link cell morphology, bioenergetics and sensory adaptation with emergent movement patterns. Exploratory analysis of the F1 trait landscape identified a varaible mixture of individual-level and cross-dependent effects, including substantial variation in growth rate between individuals and strong effects of different crosses on morphology and motility. Experimental diatom breeding may offer a unique strategy to study ocean protist evolution and phenotypic diversification and could complement other biotechnological innovations to enhance cultivation yields and crop resilience in mass cultivation.

The potential for hybridization to act as a driver of genetic diversity and adaptation in harmful algal bloom-forming species has received scarce attention, despite growing recognition of its occurrence in diverse protist and algal lineages. Prymnesium parvum s.l. (Haptophyta), is a cryptic species complex whose members form ecosystem-disruptive toxic algal blooms around the world. A prior genome analysis showed that UTEX2797, a widely used laboratory strain, originated via hybridization between clade A1 and clade A2 of this species complex. To assess the extent of A1xA2 hybridization in P. parvum s.l., we screened the genomes of 28 strains and identified 16 additional A1xA2 hybrid strains isolated from inland Texas or the eastern United States between 2001 and 2020. Chloroplast haplotypes indicated that hybridization between A1 and A2 may have occurred multiple times, and hybrids with different chloroplast haplotypes have been co-isolated from blooms in Texas in 2013 and 2020. Additionally, strain NIES1812 from Okinawa, Japan was sufficiently divergent from A1 and A2 to warrant designation as a separate clade, which we name A3. These results provide evidence for a facultative sexual life cycle in P. parvum s.l. and expand our understanding of the extensive cryptic genetic diversity present in the species complex. The frequent isolation of hybrid strains from North American blooms suggests that hybridization is common and represents a significant source of adaptive potential in these economically and ecologically damaging organisms.

The mechanisms driving host-symbiont associations across space and time in contemporary mutualisms can give insight into the capacity for symbiotic organisms to respond to environmental change. High specificity between partners can increase cooperation and facilitate efficient holobiont selection, whereas low specificity could reduce host benefit, but facilitate adaptive associations across heterogeneous environments. The present study explores specificity in natural populations of a cnidarian-algal model, Exaiptasia diaphana, across a latitudinal gradient to understand the genetic and environmental effects driving host-symbiont associations, and their relation to heritable and/or environmental symbiont acquisition. We found that symbiotic associations were extremely flexible in E. diaphana, regardless of transmission mode. E. diaphana were capable of associating with diverse symbiont communities across genetically identical hosts seeded with vertically transmitted symbionts, as well as across highly connected host populations which acquire symbionts horizontally. Host population connectivity was complex and unrelated to geographic distance, whereas symbiont community composition tracked the thermal gradient, potentially due to context dependent biotic interactions. These results indicate that in a flexible symbiosis, symbiont communities are environmentally-determined, suggesting the future of this symbiosis will likely depend on climate adaptation of symbionts.

Compared with bacterial and archaeal extremophiles, single-celled eukaryotes living in extreme habitats are understudied and underrepresented in genomic databases. An exception is the obligately halophilic stramenopile Halocafeteria seosinensis strain EHF34. A transcriptome-focused analysis of this extremophilic protists revealed the importance of organic osmolyte regulation and transport in its adaptation to hypersaline environments. However, genomic resources for H. seosinensis are currently limited to a highly fragmented assembly generated by short-read sequencing, which has hindered further investigation of the genome biology and evolution of this fascinating organism. Here, we used long-read Oxford Nanopore sequencing to generate a highly contiguous, chromosome-scale genome assembly for H. seosinensis. The assembly is 38.8 megabase pairs (Mbp) in size and contains 60 nuclear contigs, making it the most contiguous genome for a member of the order Bicosoecida. Approximately 19% of the genome is comprised of transposable elements. Of the 11,684 predicted protein-coding genes, many appear to be associated with DNA mobility-related functions, and several may be linked to adaptation to a hypersaline environment. Analysis of the H. seosinensis long-read genome assembly presented herein will facilitate our understanding of the ways in which protists have adapted to extreme environments. SignificanceHalocafeteria seosinensis is an extremophilic protist adapted to hypersaline environments. Previous analyses of a transcriptome and short-read draft genome assembly for this organism provided insights into the molecular mechanisms underlying osmotic regulation, which facilitate its adaptation to high-salt conditions. However, the lack of contiguity and quality of the draft assembly prevented the characterization of complex genomic regions, including transposable elements and viral insertions, as well as genomic comparisons with related species. Here we present a highly contiguous, chromosome-scale genome assembly for H. seosinensis that enables accurate gene prediction, detailed analysis of repeat content, and comparative genomic analysis. This long-read genome assembly will serve as a valuable resource for studying one of the few tractable halophilic protists sequenced to date.

The maintenance of endosymbiosis in cnidarians depends on the tight regulation of host immunity, cell cycle, and nutrient exchange, yet how these processes are impacted by interacting environmental stressors remains largely unknown. To address this, we employed physiological metrics, gene expression analysis, microbiome characterization, imaging (NF-{kappa}B localization, endoplasmic reticulum ultrastructure, EdU labeling), and stable isotope tracing in the model sea anemone Exaiptasia diaphana to examine the effects of heat and nitrate on these regulatory processes, individually and in combination. Heat treatment led to NF-{kappa}B activation, proteostatic stress, suppression of nutrient exchange, decreased cell-cycle progression, and microbiome restructuring, with all effects more pronounced in symbiotic than aposymbiotic anemones. In symbiotic anemones, nitrate partially offset these heat-induced responses through sustained carbon translocation, suggesting that the presence of symbionts, in conjunction with elevated nitrate, can temporarily buffer host thermal stress. However, prolonged combined exposure resulted in holobiont failure. These findings reveal that while nitrate enrichment can transiently delay the onset of bleaching, it does not preserve the regulatory networks required for symbiotic stability -- underscoring the vulnerability of cnidarian holobionts to the compounding effects of warming and nitrate pollution.

The efficiency of marine diatom carbon concentrating mechanisms (CCMs) play a critical role in photosynthesis and enable cells to maintain rapid growth rates under a variety of environmental conditions. To assess the plasticity of the model diatom P. tricornutums CCMs, specifically carbonic anhydrase (CA) enzymes and bicarbonate transporters, we measured growth response, bulk CA activity, and corresponding shifts in the proteome under a range of Zn and pCO2 conditions in culture. CA activity increased with Zn availability and decreased with pCO2. A positive growth effect was observed due to Zn addition and increasing pCO2 from 200 to 400 ppm, however growth rate decreased as pCO2 further increased to 750 ppm. Across the six treatments, the protein abundance of ISIP2A, which functions to bring Fe into the cell via a FeCO3 complex and is used as a biomarker for Fe stress, demonstrated an inverse relationship with [CO32-], consistent with its role as a phytotransferrin. Under conditions of Zn limitation ([Zn2+] = 0.3 pM), the cell appeared to allocate this metal away from CA, instead relying on a Mn-CA with a 100-fold lower intrinsic activity than that of the primary Zn-CA, as calculated using paired abundance-activity measurements. We further observed a continued increase in bicarbonate transport protein abundance after CA activity plateaued at 1.2x10-6 (reactions sec- 1cell-1), suggesting any deficit in DIC required to maintain high growth rates is accomplished through HCO3- uptake. We hypothesize that bicarbonate uptake and CO2 diffusion operate in tandem via CA enzymatic activity to supply adequate CO2 for photosynthesis.

Primary production in the Southern Ocean (SO) plays a critical role in regulating the global carbon cycle, yet the physiological mechanisms governing phytoplankton responses to iron (Fe) limitation and variable light remain poorly constrained. Using a custom made incubation system that simulated natural diel solar variability, we examined the interactive effects of Fe availability, light intensity, and photoperiod (continuous vs. variable) on three ecologically important SO phytoplankton: Fragilariopsis cylindrus, Phaeocystis antarctica, and Thalassiosira antarctica. Physiological, photophysiological, and proteomic measurements revealed that Fe availability was the dominant factor regulating growth, carbon production, photosynthetic performance and protein expression across all species. Distinct acclimation strategies emerged: F. cylindrus exhibited marked trade-offs between productivity and photoprotection under Fe stress, consistent with adaptation to stable, low-light, Fe-poor environments; P. antarctica maintained growth by flexibly modulating photoprotective and photosynthetic capacity, reflecting high plasticity suited to dynamic, open-ocean conditions; and T. antarctica expressed a balanced strategy, sustaining productivity and photoprotection simultaneously, characteristic of coastal bloom formers with higher Fe demand. Dynamic light regimes produced smaller, species-specific effects, influencing chlorophyll content and carbon storage primarily in T. antarctica. Correlation and z-score analyses demonstrated that Fe-rich photosynthetic proteins co-varied with biomass production, whereas photoprotective traits clustered independently, underscoring divergent energy-allocation strategies. Together, these results reveal how SO phytoplankton partition resources between productivity and photoprotection under shifting Fe-light regimes, providing mechanistic insight into their ecological niches.

Green algae must continuously balance resource availability to maintain photosynthetic performance. The O2:CO2 ratio is a key determinant of their metabolic mode. Under hyperoxia or low CO2, many algae induce a carbon concentrating mechanism (CCM). In the model green alga Chlamydomonas reinhardtii, the CCM relies on a pyrenoid, a specialized microcompartment that elevates CO2 around rubisco. While ambient CO2 acclimation is well-studied, responses to hyperoxia remain poorly understood, despite its frequent occurrence in nature under high light. Using controlled bioreactors, we exposed two diverse Chlamydomonas ecotypes, CC1009 and CC2343, to 95% oxygen to analyze time-dependent, genome-wide transcriptomic and phenotypic changes. Both ecotypes induced CCM genes, but they exhibited distinct molecular and physiological phenotypes. The tolerant ecotype (CC1009) successfully adapted, developing a functional CCM with a structured starch sheath. Conversely, the sensitive ecotype (CC2343) suffered growth arrest and formed malformed pyrenoids. Transcriptomics revealed that CC1009 initiated a rapid initial response, upregulating chloroplast proteostasis and downregulating nucleotide metabolism. CC2343 showed a massive, delayed transcriptional response, downregulating genes coding for photosystems and tetrapyrrole biosynthesis. This unbiased transcriptomic approach identifies key candidate genes driving algal acclimation to hyperoxic stress in natural, high-light environments.

Freshwater salinization is of increasing concern for integrity and functioning of freshwater habitats worldwide. Experiments so far often have studied drastic salt additions, while gradient designs have been performed less commonly. We tested the effect of freshwater salinization in a mesocosm exposing the plankton community of the oligotrophic Lake Lunz, Austria, to a four-fold salinization gradient (control, 0.2, 1, a 5 ppt salt). Salinity was manipulated in a factorial design with enrichment, with 10 g L-1 and 30 g L-1 phosphorus, resulting in 8 treatments with 3 replicates each. We followed the effects of salinization on diversity, community composition and resource use over 36 days. Community composition was assessed by amplicon sequencing, Diversity loss and community turnover followed upon salt addition. All levels of salinization caused pronounced changes in community composition, with 5 ppt causing the most drastic changes. Salinization caused trophic downgrading by kicking out especially protistan consumers and rotifers, while some green algae and chrysophytes were especially tolerant, resulting in reduced phylogenetic and functional diversity with increasing salinization. In line with reduced top down control, salinization affected temporal variability in chlorophyll-a (chl-a) and resource use (RUE), with higher salinity causing more extreme fluctuations in chl-a and RUE. Enrichment overall aggravated salinization, enhancing temporal turnover and temporal fluctuations in resource use.

The impacts of a changing abiotic environment on fitness and performance arise not only from low tolerance to new environmental conditions, but also from changes in the abundance and interaction intensity with other species. The strength of the interaction may itself depend on how well each species performs across environments, but there is a dearth of studies investigating how intrinsic fitness and interaction intensity covary across an abiotic environmental gradient. We addressed this question in a hypersaline consumer-resource system: the microalga Dunaliella spp. grazed by the brine shrimp Artemia franciscana. We exposed four Dunaliella strains to a range of salinities above seawater, with or without brine shrimps, and tracked their population sizes over time and the survival of their predators, to estimate basic parameters of a Lotka-Volterra model. We found that the intrinsic growth rate of algae, the survival rate of predators, and the per-capita predation rate, all varied with salinity and algal strain. Significant interactions between strain and salinity further revealed that these ecological responses to salinity are evolvable. Together with correlations between demographic parameters across salinity, this suggests that predation may influence the evolution of salinity tolerance curves, blurring the line between the fundamental and realized niches.

The symbiosis between algae and animals represents a relatively recent evolutionary innovation, and exemplified by species in the Cnidaria phylum. Cnidaria can harbor algae within a modified cellular organelle called the symbiosome in a process called endosymbiosis. This animal-algal symbiosis can be facultative or obligate. Algae acquisition occurs either by horizontal transmission, where free-swimming planula gain algae through feeding, or by algae deposition into the developing oocyte in vertical transmission[1-5]. Most studies focus on anthozoans that perform obligate endosymbiosis and transmit algae horizontally. How facultative endosymbiosis in combination with vertical algal transmission impacts the evolutionary adaptation between host and symbiont remains unclear. By studying Hydra viridissima, which performs facultative endosymbiosis and transmits its Chlorella algae vertically, we define different cell types and identify the endoderm cell lineage that gives rise to three major cell types hosting algae. Compared to obligate endosymbiosis[6, 7], Hydra viridissima algal host cells exhibit distinct features, including algal uptake, elevated oxidative phosphorylation and redox activities, and express genes that can provide ammonium for their algal symbionts. We further show when and how the developing Hydra oocytes may take up algae and where oocytes may obtain lipids. Since Hydra is amenable to genetic manipulations, our findings should enable mechanistic studies of how facultative endosymbiosis and vertical transmission evolve and adapt in a changing climate.

Near-shore coral reefs in southern Guam (Mariana Islands) experience severe sedimentation, in particular during the wet season when rainfall and erosion are high. We sampled fragments of the reef-forming coral Porites lobata from opposite ends of a sedimentation gradient in Fouha Bay, southern Guam, during dry and wet seasons. Using DNA metabarcoding, we characterized the diversity and composition of P. lobata-associated Symbiodiniaceae and bacterial microbiome communities. As in many species of Porites, Symbiodiniaceae communities of P. lobata were dominated by variants of Cladocopium C15 with sites showing differences in Symbiodiniaceae communities attributable to variation in these Cladocopium C15 variants. Bacterial microbiomes of P. lobata were dominated by Endozoicomonadaceae, a family of putative coral bacterial endosymbionts involved in nutrient cycling. Site and seasonal differences in bacterial diversity and community composition were apparent. In close proximity to the mouth of the river draining into Fouha Bay, bacterial diversity was highest during the wet season when sedimentation is generally severe. Microbiome reorganization in response to sedimentation may explain this result, but we also found overrepresentation of bacteria associated with terrestrial origin close to the river mouth and/or during the wet season. Together these patterns highlight that coral Symbiodiniaceae and bacterial communities are both spatially and temporally structured in this disturbed system. IMPORTANCEThis study provides a time series dataset of coral-associated microorganisms, including dinoflagellate algae and bacteria, from a tropical bay impacted by sedimentation that results from upstream erosion of disturbed soils. Characterizing temporal patterns of coral-associated microbes provides insights into the dynamic nature of these communities. While microbiome variability across sites and seasons may be a result of acclimatization to different environmental conditions, we identified bacterial groups of putative terrestrial origin in sampled coral microbiomes that may have been exported from eroded soils to the near-shore reef. Considering that disturbed soils act as hotspots for the proliferation of potentially harmful substances, such as antimicrobial resistance genes, understanding microbial community connections at the marine-freshwater-terrestrial interface is an important step toward evaluating environmental impacts across connected ecosystems from ridge to reef.

The non-indigenous brown macroalga Rugulopteryx okamurae has emerged as one of the most aggressive marine invaders in European waters, deeply altering benthic communities and causing severe socioeconomic impacts. While its expansion has been extensively documented along the southern Iberian Peninsula, understanding the dynamics of its northward range expansion along the Spanish Mediterranean coast remains critical for coastal management. This study documents the first formal record of R. okamurae in Calpe (Alicante), representing its current northernmost distribution limit within the Comunitat Valenciana. Sampling was conducted through an initial opportunistic scuba diving observation along the surrounding waters of the Penyal dIfac Natural Park, followed by targeted underwater surveys and an ad hoc inspection of commercial bottom-trawling nets drying at the port of Calpe during June 2026. Morphological and anatomical identification was confirmed through cross-sections of the thallus under optical microscopy, revealing the presence of both the thick and intermediate morphotypes of the species. The collected specimens were found either entangled within a native photophilic algal canopy in shallow waters or recovered from deeper offshore fishing grounds. Given the absence of records in the area during 2023-2025 surveys, these findings suggest either a very recent front-wave colonization event or a contribution from nearby, yet undetected, established patches, driven by secondary local dispersal mechanisms such as drifting fragments and explicitly highlighting commercial fishing activities as an active vector. Furthermore, considering that the species was recorded within a marine protected area and deeper environments, these results highlight a potential ecological threat to local benthic ecosystems, emphasizing the urgent need f or competent authorities to implement spatiotemporal monitoring and public awareness campaigns to prevent the definitive establishment of this invader.

Environmental DNA (eDNA) has the potential to be a powerful tool in blue carbon science for characterizing and quantifying the contribution of marine macrophytes; but its complex, dynamic relationship with bulk biomass is poorly understood. Here, we used eDNA to examine the degradation dynamics of sugar kelp (Saccharina latissima) in muddy, anaerobic marine sediment. This involved three 16-week incubations; with additions of lyophilized sugar kelp alone, a mix of lyophilized marine macrophytes including sugar kelp, and sugar kelp holdfasts buried in sediment. We used species-specific digital polymerase chain reaction assays for mitochondrial, chloroplast and nuclear markers, and metabarcoding for the 16S and 18S ribosomal RNA genes. In the former two incubations, all sugar kelp eDNA markers showed rapid log exponential declines (up to 98-99%) to asymptotes greater than the unamended controls, even as part of a more complex mix of macrophytes. In contrast, for the buried kelp holdfasts, sugar kelp eDNA increased to an asymptote (by up to [~]15X), which may be reflective of the different nature of added biomass. Overall, we demonstrate substantial preservation of environmental DNA and total organic carbon under anaerobic conditions, and the potential to use environmental DNA to quantify biomass in a blue carbon context.

Reef restoration practitioners aim to preserve coral genetic diversity by protecting reefs and cultivating diverse genotypes in coral nurseries. However, cryptic genetic lineages in most corals complicate restoration strategies, as the role of between-lineage genetic divergence remains unclear regarding adaptation. In Montastraea cavernosa, researchers have identified cryptic lineages, some strongly segregated by depth. We conducted a ten-week reciprocal transplantation experiment using two cryptic lineages restricted to shallow water (<10m depth), with one lineage more common on nearshore reefs and the other on offshore reefs. We aimed to quantify lineage-specific responses to the environment that explain the genetic and ecological divergence between the two lineages. Surprisingly, the strongest response was not lineage-specific. Instead, both lineages exhibited strong and similar changes in growth and metabolomic profiles, depending on the transplantation habitat. These results suggest that cryptic lineages employ similar mechanisms of adaptation and acclimatization to environmental challenges, despite their genetic distinction.

Phytoplankton encounter a range of light and nutrient conditions in nature and must adjust their internal carbon and nitrogen allocations to grow across different resource environments. Current phytoplankton carbon budget models treat respiration simply as a carbon loss. In reality, respiration is a critical cellular process that produces energy for nutrient uptake and cellular maintenance. Drawing on empirical evidence, we developed an eco-physiological model that incorporates a more realistic role of respiration. In our model, photosynthetic carbon is partitioned into: (i) the Pentose Phosphate Pathway (PPP) for assimilation and (ii) respiration for energy production that is then used in nutrient uptake. Stored nitrogen is partitioned between three pools: cellular structure, photosynthesis and nutrient uptake. Using an optimality-based approach, we identify strategies that maximize either exponential growth rate or competitive ability. We find that optimal internal allocations follow a growth-maintenance tradeoff, favoring population growth through carbon acquisition in nitrogen-replete conditions and population maintenance through nitrogen acquisition in nitrogen-limited conditions. The optimal allocations match empirically observed shifts in carbon partitioning at different dilution rates. Our model also generates an interactive growth response surface with an asymmetry, where light is the dominant limiting factor at low light intensities and co-limitation by light and nitrogen only occurs at high light levels. Furthermore, the model recovers the widely accepted Droop function for growth vs nitrogen quota and predicts a hyperbolic decline in growth vs energy quotas. Through a simple growth-maintenance tradeoff, our model provides a mechanistic foundation for predicting phytoplankton productivity in biogeochemical models.